Pair distribution function in a two-dimensional electron gas 



m 
o 
o 



Juana Moreno and D. C. Marinescu 
Department of Physics, Clemson University, Clemson, SC 29634 
(Dated: February 1, 2008) 

We calculate the pair distribution function, g(r), in a two-dimensional electron gas and derive a 
simple analytical expression for its value at the origin as a function of r 3 . Our approach is based 
on solving the Schrodinger equation for the two-electron wave function in an appropriate effective 
potential, leading to results that are in good agreement with Quantum Monte Carlo data and with 
the most recent numerical calculations of gr(0). [C. Bulutay and B. Tanatar, Phys. Rev. B 65, 
195116 (2002)] We also show that the spin-up spin-down correlation function at the origin, gn(0), 
is mainly independent of the degree of spin polarization of the electronic system. 

PACS numbers: 71.10.-w,72.25.-b,71.45.Gm 
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I. INTRODUCTION 

There has recently been a growth of interest in study- 
ing the pair distribution function, g(r), in electron gas 
models^*2i&4 caused mainly by its relevance in non-local 
density functional theories The zero inter-electronic 
distance value, g(r = 0), also appears in the large 
wavevector and the high frequency limits of the electronic 
charge and spin response functions.^* The importance 
of g(r) lies in its connection with the electronic exchange 
and correlation of the electron gas model. Moreover, the- 
oretical calculations of the pair distribution function can 
be directly compared with material properties since g(r) 
is the Fourier transform of the static structure factor. 

The pair-distribution function is the probability of 
finding a pair of electrons at a distance r from each other. 
Therefore, the average number of electrons in a spheri- 
cal shell centered on an given electron is ng{r)Vlo''' D ~ l dr, 
where ^Ijyr 0-1 ^ is the volume of the D-dimensional shell 
and n = N/V is the uniform electron density. At large 
distances, g(r) approaches 1, whereas near the origin, 
where the electron charge is depleted, it is small on ac- 
count of the Pauli exclusion principle and the exchange 
and correlation effects associated with the Coulomb in- 
teraction. 

The subject of this paper is an analysis of the pair cor- 
relation function dependence on the inter-electronic dis- 
tance and electron density in a two dimensional, interact- 
ing, spin polarized electron system. Calculations of g(r) 
in the two-dimensional paramagnetic electron gas have 
been reported by FreemaniS and Nagano et al. 11 within 
the ladder approximation, by Tanatar and Ceperlej^ us- 
ing the diffusion quantum Monte Carlo method (QMC) 
and, more recently, by Bulutay and Tanatar— using the 
hypernetted-chain approximation (CHNC). Moreover, an 
analytical expression of g(0) has previously been derived 
by Polini et al£ 

In order to calculate g(r), we follow the approach de- 
veloped in three dimensional systems by Overhauser— 
and further refined in Ref. and 13. This method is 
based on the relation between g(r) and the two-electron 
scattering problem in an appropriately chosen effective 
potential which will be discussed in detail in the follow- 



ing sections. In addition to obtaining the variation of the 
pair distribution function as a function of the coupling- 
strength, r s , and of the spin polarization, we also derive 
an analytic expression for <?(0): 
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This expression is found to agree very well with the re- 
sults of the most recent numerical calculationspi^iii We 
also compare our results with the expression derived in 
Ref. H 



II. EFFECTIVE MODEL FOR THE PAIR 
DISTRIBUTION FUNCTION 

A spin-polarized electron gas is characterized, in equi- 
librium, by two parameters: the electronic density, or its 
equivalent r s , and the polarization, ( = (n-f — nj_)/(7ij + 
n I ) , where and n± are the spin-up and -down electron 
densities. For this system, the pair distribution function 
is given by 

g(p) = J[(i+C) 2 3tt(p)+2(i-C 2 )3ti(p) + (i-0 2 3 U (p)], 

where g aa ' are the spin resolved pair distribution func- 
tions. Following Overhauser fiSi g aa > can be related with 
the two-electron wave functions as: 3 

9U(P) = 75<l*singlct(p)| 2 )p u (fc) + 75(l*triplot(p)| 2 )p n (fe),(3) 



9U(P) = (l*triplct(p)| 2 )p TT (fc), 
911(P) = (l*triplct(p)| 2 )p u (fc), 



(4) 
(5) 



where * s ingict(p) and * triplet (/>) are, respectively, the 
two-electron wave function for the singlet and triplet 



, (j.) denotes the average over the prob- 



states and (. 

ability of finding two electrons with relative momentum 
k and spins a and a . 3 
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The wave function of an electron pair, \&(p), verifies 
an effective Schrodinger equation: 



<9 2 * f9* f <9 2 * 



2m* V dp 2 p dp p 2 



V(p)^ = E^, (6) 



where V(p) is the effective potential, m* = m/2 is 
the reduced mass and E is the energy of the elec- 
tron pair, which is approximated by h 2 k 2 / (2m*). Since 
the solution to this equation can be written as * = 
cos (m</>)\E' m ( j o), the spin resolved pair distribution 

m 

functions become 

oo 

9U(P) = (l*o(p)| 2 ) m(fe) + 2 J2 (\* m (p)\ 2 ) Pn (*), (7) 

m— 1 

5TT(P)=4 (\^n(p)\ 2 ) Pn{k) . (8) 

m odd 

Overhauser's method relies on the appropriate selec- 
tion of an effective potential capturing the short range 
correlation effects of the Coulomb interaction. In three 
dimensions, Overhauser chose the electrical potential cre- 
ated by an electron and a neutralizing sphere of uniform 
charge with radius r s surrounding it4^*i^ The effective 
potential is expected to mimic the true one when the rel- 
ative distance between electrons verifies r < r s . When 
r > ?- s the potential vanishes and is not expected to be 
close to the true potential felt by an electron moving 
in a uniform electron gas. This approach is equivalent 
to assume that the probability of finding three electrons 
in a sphere of radius r s is exactly zero£ Numerical es- 
timates of this probability for a three-dimensional inter- 
acting electron gas^ have shown that is indeed small and 
we expect the same result holds in two dimensions. 

Following this procedure, in two dimensions, we might 
approximate the screened Coulomb potential by the po- 
tential of an electron surrounded by a circle of radius 
r s uniformly filled with screening charge density ne — 
e/(irr 2 ). For convenience, we introduce dimcnsionlcss 
variables, x = p/r s and V(x) = V(p)/(e 2 /r s ), where 
r s is measured in units of Bohr radius (as — h 2 /me 2 ), 



V(x) = -- -E(x), 

X 7T 
f 4 

V(x) = x E 

X 7T 



(l) 



1 



x € 1 



AW 1 

x* / \x 



(9) 
x > 1 



where K(x) and E(x) are, respectively, the complete el- 
liptic integral of first and second kind. The screened 
potential of a uniformly charged disk of radius r s with 
an electron at its center does not vanish, but it has an 

attractive long-range tail, V(x — ► oo) — > -. Since 

we are interested on obtaining an analytical expression 
for g(0), a further simplification of the effective potential 
is needed. Since Overhauser's effective potential is not 
reliable outside the disk of radius r s , the most reason- 
able simplification is to make it zero outside this disk. 



0.4 



0.2 



3 o 
> 

-0.2 



-0.4 



\ 


v 6 „ 

— v 6H -v 

I Polini et al. 

\ 
\ 
\ 
i 
\ 
■■\ 


\ 


/ — " 


\ \ 
\ \ 
\ x 


/ ^ — 



0.5 



p/r s 



1.5 



FIG. I: Normalized effective potential, V(p)/ (e 2 /r s ), as func- 
tion of p/r s . The potential from Eq. (solid line) and our 
choice of effective potential, Eq. (1101 . (dashed line) are dis- 
played. The effective potential used by Polini et all— is also 
shown (dotted line). 



To avoid a discontinuity in the effective potential and 
considering that V(x) is arbitrary to the extent that a 
constant can be added to it, we subtract from the po- 
tential in the region where x ^ 1 its value at x = 1, 
Vq = 1 — (4/7r). Thus, our effective potential is: 



f 4 4 
V(x) = E(x) + - 

X 7T 7T 

V{x) =0, x ^ f . 



f, X < 1 



(10) 



Figure ^ displays the initial effective potential from 
Eq. © together with our election of effective potential, 
Eq. Ijf and Polini et al. choice 4 , which was based on 
a previous variational calculation^ The main difference 
between the effective potential used in Ref. and ours is 
that the former one has a discontinuity at p = \pKr s j2 
while ours is always continuous. 

Using Eq. (|1(JI) for the electronic potential, the 
Schrodinger equation becomes: 



d 2 V T 
dx 2 



1 d^ r 



x dx 
d 2 *„ 



7TI \ 

V 2 -— J*m(aO = , a? > 1, 



f d^, 



-r.„ - 



dx 2 
4 



E(x 



x dx 

7T 



* m (x) = o, x s; i,(ii) 



where the relative momentum is also renormalized, q — 
kr s . The general solution for x f is given by ty m (x) = 
Jm(qx) + B m (q,r s )N m (qx), where J m is the Bessel func- 
tion of order m and N m is the corresponding Neumann's 
function. The coefficient B m (q,r s ) can be written as 
B m (q,r s ) = cot(S m (q,r s )), where S m (q,r s ) is the wave 
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rent relation between the coefficients: 



FIG. 2: Pair distribution function of the unpolarized electron 
gas at r s = I as a function of rkF. Our approximation (solid 
line) is compared with the Quantum Monte Carlo data of 
Ref. ^3 (crosses) . 
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FIG. 3: Pair distribution function of the unpolarized electron 
gas at r a = 5 as a function of rkF- Our approximation (solid 
line) is compared with the Quantum Monte Carlo data of 
Ref. ^3 (crosses) . 



function phase shift due to the presence of the scattering 
potential^ 

To find the solution inside the circle of radius unity 
we make a Taylor expansion of the pair wave function: 

oo 

^ m (x) = OL m ^ l x n . We arrive to the following recur- 



(n 2 - m 2 )a m ^ n = r s |a m:n _i + ( 3^a mj „_ 2 

(12) 



where A(r) = 2 



(2r- 1)!! 



2 r r\ 



1 



2r- 1 



As a consequence 



of this recurrent relation, every a m) „ is proportional to 
a mim and a function of r s and q, a m ^ n = a m . m F n (r s , q). 

In order to solve Eq. we match ^ m (x) and its 
derivative at x = 1: 

G m (r s ,q) = J m (q) + B m (q, r s )N m (q) (13) 
F m (r s ,q) = qj m (q) + B m (q, r s )qN m (q) (14) 

OO 

where G m (r s ,q) = ^ F n (r s ,q) and F m (r s ,q) = 

n—m 

oc 

nF n (r s , q). For a given momentum transfer and cou- 

n—m 

pling strength the parameters a m . m and S m (q, r s ) become 

/ \ J m {q) +cot(S m (q,r s ))N m (q) 
a m . m (q,r s ) = — , and (15) 

^m^s, q) 



cot(S m {q,r s )) 



(r s ,q)J m (q) - G m (r s ,q)qJ m (q) 



G m (r s ,q)qN m {q) - F m (r Sl q)N m {q) 

The pair wave functions ^ m {x) are computed for any 
value of q and r s using Eqs. (fl5|l and iflfUl and the spin 
resolved pair distribution functions are calculated using 
Eqs. J7J and © and an appropriate choice for the dis- 
tribution of the relative momentum of an electron pair. 
For simplicity, we use the probability distribution of a 
free Fermi gas. For the unpolarized electron system, the 
probability of a pair with momentum q is independent 
of the spin orientation and proportional to the overlap 
between two circles of radius k F displaced by 2q^ 



Paa 1 (<?) 



16<? 
nk 2 F 



arccos 



(17) 

Figures [21 and [3] display our results for the pair distri- 
bution function of the unpolarized system at r s = 1 and 
r s = 5, respectively. We have used seven (m max — 7) 
partial waves and up to n = 50 terms for the expansion 
of \E' m (a;) in the internal disk. Our results at moderate 
coupling strengths agree quite well with the QMC data. 12 
However, at larger values of r s our method is unable to 
reproduce the strong quantum oscillations of the numer- 
ical results. This discrepancy is expected and shared 
by previous calculations using a self-consistent Hartree 
schemeJ. With decreasing dimensionality the role of ex- 
change and correlations becomes more important and a 
screened Coulomb potential is insufficient to completely 
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capture this physics. The introduction of self-consistent 
spin-dependent effective potentials have probed able to 
reproduce more closely the numerical results in this range 
of densities^ 



III. PAIR DISTRIBUTION FUNCTION AT THE 
ORIGIN 

At zero distance, g-\-\{p = 0) vanishes on account of the 
Pauli exclusion principle, while g-\i(p = 0) is determined 
by the m = component of the two-body wave function, 



ffu(p = 0) = (\^o(p- 



)pu(*0 = (l«o,o| 2 ) Pn (fe)- (18) 



Since the distribution of the relative momentum of an 
electron pair is a smooth function, a good estimate of 
(l a o,o| 2 )p n (fc) is obtained by making an expansion around 
the momentum where the distribution reach its maxi- 
mum as: 



1 



«o,o 



G(r s 



(19) 



where the momentum dependence of G have been 
dropped. Using the recurrent relation (|12|) . we obtain 
a series expansion of G(r s ): 



oc 

G(r s ) = l + r s { ai +J2 



A(m) 



OO 

:{a 2 + £x(m)[- 



^ (2m + 2) 
1 



m— 1 

4/ (2m + 4) : 



2m + 3) 2 



1 



1 



(2m + 2f 



(2m 



i 



(2n + 2) 2 (2n + 2m + 4) 5 



}+0{rl) 



- (1 + 0.6032 r s + 0.07263 r 2 s ) 
where a\ 



(20) 



1 1 5 11/13 

— I and a9 = 1 — 

4 7T - 9tt 6 4 \tt 4, 

We can obtain g(0) at any order in the expansion 
on the parameter r s since that g(p = 0) = — 5tx(0) = 

— — — — . To first order in the expansion of G(r s ) the 
2 G(r s ) 

pair distribution is: 



/(p = Q) 



2 [1 +0.6032 



(21) 



To second order we recover Eq. Our approximation 
procedure also allows us to sum all the orders as: 



g{p = o) 




(22) 



Finally, we also calculate g(0) performing the average 
over the distribution of relative momenta: 



1 



3(0) = 7T5Tl(0) 



p(q)\a ,o(q,r s )\ 2 dq (23) 



where o>ofl(q,r s ) is given by Eq. (|15f) and l|16l) and p{q) 
by Eq. (|T7|l. 

Fig. 01 displays the pair distribution function of a two- 
dimensional unpolarized electron gas at the origin as a 
function of the coupling-strength r s . Results of the first 
order in the analytic expansion on r s , Eq. I|21l) , the second 
order, Eq. J[J, and the infinite order solution, Eq. I|22|) . 
together with the momentum average results, Eq. 12311 . 
are displayed. In addition, the results of the numeri- 
cal calculation of Bulutay and Tanatar— and the recent 
estimate by Polini et alX are also included for compar- 
ison. Several conclusions can be gathered. By adding 
additional terms in the analytical expansion on r s we are 
able to closely approach the numerical results^ in the 
low-density regime. Note that Eq. is already a reli- 
able analytical expression for g(0). Fig.0]also shows that 
the results of the momentum average approach, Eq. 123|l . 
are slightly below Bulutay and Tanatar results for small 
r s , but become even closer to the numerical curve in the 
low-density regime. In this regime, the analytical ex- 
pression obtained on Ref. displays much larger values 
of g(0) than the available numerical dataiiSiH 

The contact value of the pair distribution function 
changes when the electron gas is polarized. The spin 
polarization directly appears on the expression for <?(0), 

1 - c 2 

g(0) = — - — <7||(0). In addition, the polarization modi- 
fies the distribution of momenta of the electron pair and, 
as consequence, the value of g-\i(0). We calculate the spin 
resolved pair distribution function as: 



3Ti(C,P = 0) 



Pd<l)\ a o,o(q,r s )\ 2 dq. 



(24) 



where p^(q) is the distribution of relative momentum in 
the polarized electron gas: 



pc(<z) = 



(for ^ q ^ fc_) 



pc(i) 

+k 2 



Sq 



irk 2 k 2 



I ( arccos(x) — x-\f\ 



F[ ( arccos(y) - yy/l - y 2 

(for fc_ q sC k+) 



(25) 



where x 
k 2 - k 2 



.9 



k 2 



k 2 
K Fl 



4q 



Aft, V 



k i and fc_ = \kpi — fc.FjJ/2, fc+ = (k 



FT 



fcpjJ/2. The Fermi momentum for the spin up (down) 
population is related with the polarization and the 
Fermi momentum of the unpolarized gas (kp) by kp j = 
kpy/1 + £ and kpi = kp^Jl — C,. 
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FIG. 4: Two-particle distribution function at the origin as a function of r s . Results of the analytical expansion to first order 
given by Eq. 12 H (dot-dashed line), second order by Eq. Q (dashed line), and infinite order [Eq. 1221 1 (dot-double-dashed 
line), and the average over the distribution of relative momentum [Eq. (I23H 1 (long-dashed line) are displayed. The numerical 
results of Bulutay and Tanata»— (solid line) and the interpolation results of Polini et alA (dotted line) are also displayed. The 
right panel is a blow-up of the large r s region. 



Our results show that g~\i((,p = 0) is largely unaf- 
fected by the degree of spin polarization. The difference 
between g^i(C = 1)0) an d its unpolarized counterpart 
5Tl(C = 0)0) is, at most, a few percents for any given 
value of r s . The absence of a significant dependence 
with the spin polarization was also found in previous 
calculations ^ Moreover, given that the momentum de- 
pendence of our results is rather weak, we do not expect 
important changes if the free Fermi momentum distribu- 
tions, Eq. 1)17(1 and 1)25(1 . are replaced by the interacting 
ones. 



IV. CONCLUSIONS 

We have calculated the pair distribution function in 
a two dimensional interacting electron gas, following an 
approach originally developed in three dimensions by 
Overhauser^ii Within this framework, the short range 
correlations of the Coulomb interaction are replaced by 
an effective potential, and the calculation of g(p) is re- 
duced to solving the corresponding two-electron scatter- 
ing problem and averaging over the probability distribu- 
tion of the momentum of the electron pair. Our results 
for g(r) at moderate coupling strengths agree well with 
the numerical data^ At larger values of r s , however, this 



approximation is unable to reproduce the strong quan- 
tum oscillations of the numerical results. 

The analytic expression for g(0) as a function of 7' s , 
Eq. derived in this context compares very favor- 
ably with the complete solution of the effective potential, 
Eq. 1(23(1 . and with recent numerical calculations»i2iii We 
believe that the discrepancy between the present results 
and the analytical expression obtained in Ref . is essen- 
tially due to the different choice of effective potential (see 
Fig. ^| . Besides, while we have used the same approach 
for all values of the electronic density and polarization, 
Polini et al. use an interpolating scheme between the 
results of a perturbative expansion at high-density and 
the Overhauser's treatment of scattering processes in the 
low-density limit. 

We have also studied the dependence of g(0) with the 
spin polarization of the electron gas. We have found that 
the spin-up spin-down correlation function, g|x(0), is ba- 
sically independent of the degree of polarization. There- 
fore, the polarization modifies g(0) only through its de- 

2 /It Tt\ 

pendence on the density, g(0) = t. — 9u(0). 

n 

Within this approach, further study of how the choice 
of the effective potential modifies the pair distribution 
function can provide valuable insight into the short range 
electronic correlations in real materials. 
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